Local Gating of an Ir(lll) Surface Resonance by Graphene Islands 
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The influence of graphene islands on the electronic structure of the Ir(lll) surface is investigated. 
Scanning tunneling spectroscopy (STS) indicates the presence of a two-dimensional electron gas with 
a binding energy of —160 meV and an effective mass of —0.18 m e underneath single- layer graphene on 
the Ir(lll) surface. Density functional calculations reveal that the STS features are predominantly 
due to a holelike surface resonance of the Ir(lll) substrate. Nanometer-sized graphene islands act 
as local gates, which shift and confine the surface resonance. 

PACS numbers: 73.20.At, 73.21. Fg, 73.22.Pr 



On a number of metal surfaces, single layers of 
graphene may be grown. Owing to their different 
graphene-metal interactions [1|, these substrates may 
modify the electronic structure of the graphene layer. 
For example, the tt and 7r* bands of graphene, which give 
rise to a Dirac cone dispersion in the pristine material, 
are significantly altered on Ni(lll) and Ru(0001) 0,0, 
where a gap in the Dirac bands opens. In contrast, from 
Ir(lll), an almost unchanged band structure of graphene 
has been reported [4] . The modification of the electronic 
structure of the substrate upon graphene adsorption has 
hardly been investigated. From calculations, quenching 
of the Ni(lll) surface state upon graphene adsorption 
was predicted [5[ . The only experiment related to this is- 
sue addressed a spatial variation of the Ru(0001) d states, 
which arises from a graphene-induced moire pattern [6|. 
This state of affairs is surprising as electronic states at 
surfaces are sensitive probes of the interaction with ad- 
sorbates. 

Here, we combine scanning tunneling spectroscopy and 
density functional theory (DFT) calculations of Ir(lll) 
covered with a single layer of graphene. Spectra of the 
differential conductance (dl/dV; I: current; V: sam- 
ple voltage) of pristine and graphene-covered Ir(lll) re- 
veal a holelike surface resonance. Intriguingly, the most 
prominent features in the spectra are due to this reso- 
nance at the Ir-graphene interface rather than to any 
states of the graphene layer. This effect is attributed 
to the selectivity of the tunneling current for states 
with small parallel momentum. Although the electronic 
structure of graphene is only weakly perturbed by the 
Ir(lll) substrate, graphene shifts the Ir surface reso- 
nance by ~ 190 meV towards the Fermi level. As a re- 
sult, graphene islands act as local gates which confine 
the surface resonance and induce characteristic standing 
wave patterns. The resonance shift and an effective mass 
m* = — 0.18m e (m e : free electron mass) determined from 
these patterns are consistent with DFT results. 

Experiments were performed with a home-built scan- 
ning tunneling microscope (STM), operated at 5.2 K in 
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FIG. 1. Main figure: Spectra of dI/dV/(I/V) recorded along 
a line from Ir onto a graphene island as indicated in Inset (b) . 
A clear shift of the surface resonance onset occurs. Contours 
of constant LDOS calculated using a scattering model (see 
text) are indicated by black and light gray lines for Ir and 
graphene, respectively. A horizontal blue (gray) line indicates 
the position of the graphene edge. Inset (a): Spectra of dl/dV 
recorded near the beginning (blue (black)) and end (red (dark 
gray)) of the line in Inset (b). The beginning of the line is 
defined by the cross at the dotted line in Inset (b). Inset (b): 
Constant-current STM image of a graphene patch on Ir(lll) 
(-220 mV, 100 pA). White dots mark positions where dl/dV 
spectra were recorded. A white cross denotes zero distance. 
The false colors used in the main figure and in Inset (b) are 
defined in the upper-right-hand corner. 



an ultrahigh vacuum. Ir(lll) surfaces were cleaned by 
cycles of Ar + bombardment and annealing. Graphene 
films were grown by exposing the sample to ~ 6 x 
10 _4 Pas of C2H4 at room temperature and subsequent 
annealing at ~ 1400 K. This procedure leads to a par- 
tial coverage of Ir(lll) by highly ordered graphene Q. 
Au tips were prepared by ex-situ cutting, in-vacuo heat- 
ing and Ar + bombardment. Spectra of dl/dV were ac- 
quired by a standard lock-in technique (modulation fre- 
quency: 9.1kHz, modulation amplitude: 10mV rms ) and 



subsequently normalized by I /V to compensate for the 
voltage-dependent transmission of the tunneling barrier 
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Figure Q] shows normalized dl/dV data, which were 
obtained along a slightly curved line [Inset (b) to Fig. [I] 
crossing a step from bare Ir onto a graphene island. 
Above the bare Ir(lll) surface, a steplike decrease in the 
dl/dV signal at V « -350 mV [Inset (a) to Fig.ffl blue 
(black) line] occurs, which is indicative of a holelike sur- 
face state or resonance [9j . Such a resonance was previ- 
ously observed with photoelectron spectroscopy [l(J [H| • 
Above the graphene layer, the steplike feature is shifted 
to V ~ — 150 mV [Inset (a) to Fig.HJ red (dark gray) 
line] . When approaching the edge of the graphene island 
from either side, the dl/dV step moves towards lower 
voltages and disappears on top of the graphene edge. 
This spatial variation of the shift can be explained by 
scattering at the graphene edge. The rather strong inter- 
action between graphene edges and the Ir(lll) substrate 
bends the graphene edges towards the metal and leads to 
the formation of chemical bonds between Ir and C atoms 
at the island edge [12] . The Ir(lll) surface resonance is 
considered as a free electron gas with a binding energy 
Eq and effective mass m* scattered from a hard- wall po- 
tential provided by the graphene edges. The spatial vari- 
ation of the local density of states (LDOS), p s , can then 
be described as 



p s (E,x) = 1 - Jo(2fe||x), 



(1) 



where Jo is the zeroth-order Bessel function, k\\ = 
^2m*(E — Eo)/h is the parallel momentum, E is the en- 
erg y, and x = is the position of the hard- wall potential 
[13Ml5|. Lines in Fig. [T] show contours of constant LDOS 
calculated for m* = — 0.18m e (see confinement analy- 
sis below), E'ojir = — 350 me V (curved black line) and 
^o,gr = — 160meV (curved white line), for Ir(lll) and 
graphene-covered Ir(lll), respectively. To match exper- 
imental data, a Gaussian broadening of 40 meV was ap- 
plied. Further, the energies Eo jgr and i?o,ir were chosen 
to yield a good fit between the calculated first maximum 
of the oscillation (embraced by contour lines) and the ex- 
perimental data. The simple hard- wall model reproduces 
the curvature and position of the LDOS maxima quite 
well and yields an energetic shift of the surface resonance 
between bare and graphene-covered Ir of AE ^190 meV. 
The energies fit well to the energies extracted from the 
single spectra; Eo ygr also matches the energy obtained by 
confinement analysis (see below). It is important to note 
that the data do not reveal any particle-hole symmet- 
ric counterpart of these confinement features above the 
Fermi level. At energies below « — 650 me V (not shown 
in Fig. [I]) variations in the LDOS with the periodicity of 
the moire pattern predominate, probably due to weak 
periodic potential modulations. 

To further support the above model, DFT calcula- 
tions of the pristine and graphene-covered Ir(lll) surface 
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FIG. 2. Calculated band structures of (a) pristine and (b) 
graphene-covered Ir(lll). Light gray lines in (a) show the 
dispersion of all states of the supercell used. The contribu- 
tion of each state to the tunneling current is indicated by the 
widths of blue (black) and red (dark gray) lines. On Ir, the 
current is essentially due to a surface resonance at the center 
of the surface Brillouin zone (r). On the graphene-covered 
surface, the resonance is shifted upwards. It still carries most 
of the current, while the Dirac cone around K is less conduct- 
ing. Green (gray) parabolas show a fit with an effective mass 
m* « -0.17m e . 



were performed [l6|. We used the projector-augmented 
plane wave method [13, EH , as implemented in the Vi- 
enna ab initio simulation package [19]. Light gray lines 
in Fig.[2^a) show the calculated band structure of an 
Ir(lll) surface. Near the Fermi energy of the Ir(lll) 
system, various bands originate from bulk bands with p 
or d character and there are several surface states around 
the bulk band gap near K and a surface resonance near 
T. This resonance is mainly derived from Ir p z orbit als 
of the first few atomic layers near the surface. In the 
energy range — leV < E < OeV probed in the STM 
experiments, the calculated dispersion of the surface res- 
onance is approximately parabolic, with an effective mass 
of m* « — 0.17m e [Fig.[2fa), green (gray) line]. At 
the graphene-covered surface [Fig.[2fb)], two additional 
bands derived from the C p z orbitals occur and form a 
Dirac cone near K in agreement with previous photoemis- 
sion [11] and DFT [4] studies. On the graphene-covered 
surface, the resonance is shifted upwards [Fig.EJb), green 
(gray) line] by an amount which depends on the distance 
between the graphene sheet and the topmost Ir layer. 
For typical spacings between 0.327 and 0.362 nm [20|, 
the calculated shift is between 100 and 200 meV, which 
is consistent with our experimental value of « 190 me V. 
To trace back the origin of the resonance shift upon 
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graphene adsorption, calculations were performed in 
which the graphene C atoms were replaced by chemically 
fully inert Ne atoms. As a result, a Ne layer shifts the Ir 
surface resonance upwards by virtually the same amount 
as the graphene layer. As Ne provides no states at the 
Fermi level which could donate or accept charge from 
the Ir surface, the upward shift of the resonance is most 
likely due to a significant Pauli repulsion. Nevertheless, 
Coulomb potential effects, e. g., via charge redistribution 
[20|, occur and cannot be disregarded in modeling the 
full electronic structure of graphene/Ir(lll). The calcu- 
lations further show a downward shift of the Dirac cone 
when the graphene sheet is pushed towards the Ir sub- 
strate. This shift cannot be explained by Pauli repulsion 
and demonstrates that Coulomb potential effects are pre- 
dominant for the energy of the graphene Dirac point. 

To determine the contributions of the various states 
to the tunneling current, the approach of Tersoff and 
Hamann [21] was used and the tip was modeled as an s or- 
bital \L). Based on an estimate of the experimental tip- 
sample distance (22| the orbital is placed 0.48 nm above 
the surface. The overlap \ {^ n ^\L)\ 2 , where \^ n ,k) is the 
wave function of a band n at wave vector fc, is indicated 
by the width of the colored bands in Fig. El On both 
surfaces, clean and graphene-covered Ir(lll), the main 
contribution to the current is due to the aforementioned 
surface resonance. The current due to the Dirac bands 
of graphene is significantly smaller. This may be under- 
stood from the parallel momenta fey of these states, which 
affect their decay into vacuum. The surface resonance is 
located around T and thus decays less rapidly than the 
Dirac cone states near K [23|. Therefore, the steps in the 
dl/dV spectra may safely be attributed to the (shifted) 
Ir(lll) surface resonance. This result is also in agree- 
ment with the absence of electron-hole symmetry from 
the experimental spectra. In recent publications 24], 25[ , 
scanning tunneling spectroscopy data from graphene on 
Ir(lll) are attributed to tunneling from graphene states. 
However, the analyses of Refs. [24| and |25| neglect the sub- 
strate electronic states at the Brillouin zone center. In 
contrast, the present results show the importance of sub- 
strate states at T which give the dominant contribution 
to the current in our STM experiments. 

In addition to scattering at their edges, graphene is- 
lands lead to confinement of the hole states. Figure[3] 
shows an STM image of a graphene island along with 
normalized dl/dV maps recorded at constant current. 
At increasingly negative sample bias [Figs.|3jb)-(d)] the 
pattern inside the island evolves from a central maximum 
over a ring to a ring with a central maximum, as expected 
for confined states with zero, one, and two nodes, respec- 
tively. For a more detailed analysis, we model the island 
by a circular quantum dot with hard walls. The eigenen- 
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FIG. 3. (a) Constant-current STM image of a graphene island 
on Ir(lll) (-220 mV, InA). A blue (light gray) circle indi- 
cates the effective island diameter of 8.3 nm used for further 
analysis, (b)-(d) Normalized dl/dV maps of the graphene is- 
land in (a), recorded at the indicated voltages. The observed 
LDOS oscillations evolve as expected for the confinement of 
an electron gas in a quantum dot. (e) Energies of confined 
states with significant LDOS at the island center (n = 0) 
measured from various graphene islands. States with / = 1 
(open blue circles) and / = 2 (filled red circles) are resolved. 
Lines show a fit according to Eq. ([2]) with Eq — — 160 me V 
and m* = — 0.18m e , which are in good agreement with val- 
ues calculated within DFT. 



ergies E n j of a confined electron gas are [26 



E n ,l = Eq 



m*d 2 ' 



(2) 



where u n j is the Ith. root of the nth-order Bessel func- 
tion J n and d is the island diameter. Figure [3fe) displays 
the energies of the first two resonances which exhibit an 
LDOS maximum at the island center (^0,1 an d ^0,2) eval- 
uated from spatially resolved d//dV spectra of 8 nearly 
circular islands with diameters between 5 and 12 nm. 
The effective island diameters d were determined from 
an inscribed circle, which touches the island boundaries 
at the midpoint of the step edge [Fig.[3ja), blue (light 
gray) line] . Energies calculated according to Eq. (|2j) with 
E = -160 meV and m* = -0.18 m e [Fig.^e), lines] 
match the experimental data very well. 

As a final test of the model spatially resolved dl/dV 
spectra from a roundish island are compared with the 
calculated LDOS. The LDOS, p(E,x), of the surface 
resonance confined to a disk is 



p(E,x) 



n,l 



E n ,i — E 
SE 



(3) 



are solutions to the Schrodinger equation as de- 
scribed by Piatt et al. [261 ] . 

V ni i(r,ip) = J n exp(in^). (4) 
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FIG. 4. Spatially resolved normalized dl/dV spectra recorded 
along a line through a graphene island (inset). The experi- 
mental distance scale was recalibrated to account for the small 
curvature of the measurement path. Light gray lines mark 
contours of constant LDOS calculated for the circular island 
which is indicated in the inset. Blue (gray) lines indicate the 
boundaries of the model island. Inset: Constant-current STM 
image of the graphene island (—220 mV, 100 pA). White dots 
indicate positions where spectra of dl/dV were recorded. A 
white ring denotes zero distance. An inscribed blue (gray) 
circle shows the effective island diameter of 11.8 nm used in 
modeling. 



SE is a Gaussian broadening reflecting a finite lifetime 
of the states. While the broadening may depend on en- 
ergy 27], the constant broadening SE = 40 meV assumed 
here is sufficient to match the experimental observations. 
Figure [4] shows a series of 29 normalized dl/dV spectra 
measured along a line across a graphene island (inset). 
Light gray contour lines show the calculated LDOS using 
the measured diameter of 11.8 nm (blue (gray) horizon- 
tal lines). The qualitative agreement of the theoretical 
and experimental data is further evidence that graphene 
islands confine the Ir(lll) resonance. 

In conclusion, van der Waals-bonded graphene on 
Ir(lll) induces a pronounced shift in the Ir(lll) sur- 
face resonance. The disappearance of the resonance at 
graphene edges indicates the covalent carbon-metal in- 
teraction, which acts as a hard-wall potential for scat- 
tering of resonance electrons. Nanometer-sized graphene 
flakes can therefore confine quasi-two-dimensional elec- 
tron gases to artificial quantum dots. 
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Note added in proof: Results concerning confined elec- 
tronic states in graphene islands on Ir(lll) [Hj have been 
published after submission of this manuscript. Ref. [28| at- 
tributes the confined states to either graphene states or 



the scattered Ir(lll) surface resonance, depending on the 
graphene island size. This is in contrast to our interpre- 
tation that the surface resonance predominates the STM 
data for any island size. 

During the refereeing process an experimental obser- 
vation of a graphene-induced shift of the Ir(lll) surface 
resonance has been reported [29[. 
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